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Gauge theories predict the existence of topological defects, such as, do- 
main walls, cosmic strings, and monopoles, in the early stages of the Universe 
[1]. Among them, the most studied are cosmic strings (CS), because they 
may provide seeds for the formation of galaxies and the large scale struc- 
ture of the Universe [2] . This scenario of the structure formation has been 
extensively studied in the past decade [3], and some of its fundamental as- 
sumptions have been challenged by some recent numerical simulations [4]. 
In particular, it was found that there is a significant amount of small-scale 
structure on scales from the horizon down to the limit of resolution of the 
simulations. As a result, the single-scaling regime, upon which the original 
scenario rests, seems to be ruled out. Most of the small-scale structure is 
due to kinks, which are formed on the long CS when they cross and inter- 
commute. These kinks emit gravitational radiation, which is very similar, in 
both frequency and amplitude, to the one from loops [5]. 

The study of the gravitational radiation of CS is of particular importance 
to the cosmological scenario, since most of the observational limits on CS 
networks in the early Universe are obtained from it, and more important, 
the multiple-scaling regime [3] crucially depends on the rate of the radiation. 
If the radiation is strong enough so that the small-scale structure can be 
smoothed out by it, a multiple-scaling is resulted, and the scenario can be 
saved. Thus, to have the scenario work, a complete understanding of the CS 



gravitational radiation is required. 

In this Letter we shall study the backreaction of the gravitational and 
particle radiation emitted by infinitely long cosmic strings (ILCS), by con- 
sidering analytical models to the Einstein field equations. In order to avoid 
some mathematical uncertainties [6], we shall consider ILCS with finite thick- 
ness. In addition, our study of particle radiation is only restricted to the 
massless one. The main reason is that for macroscopically large CS, only 
can the emission of massless particles be important, and the ones of particles 
with masses are usually suppressed [1]. 

It should be noted that the gravitational interaction of ILCS with parallel 
moving gravitational and electromagnetic waves has been studied in [7, 8]. In 
particular, Letelier [8] showed that the interaction sometimes could produce 
intermediate (or non-scalar) singularities. However, since there is no momen- 
tum transported away from the CS, his model as well as the one studied in 
[7] cannot be considered as radiation of CS. The gravitational radiation of 
ILCS was studied in [9] by weak- field approximation (WFA). Again, because 
of the use of the WFA, the backreaction of the radiation field was not fully 
taken into account. Due to the non-linearity of the Einstein field equations, 
however, it is expected that it might play a significant role in the evolution 
of CS. In fact, in this Letter we shall show that it is this backreaction that 
always makes the spacetime of CS singular on the symmetry axis. 



To begin with, let us first consider the spacetime with the metric [10] 

where t, r, z, and (p, as cyhndrical coordinates, have their usual physical 
meaning, with — oo < t,z < +oo, < r < H-oo, and < ip < 2tt. And 
6(77) is given by b{ri) = [B{ri) — B{rio)]H{ri — rjo), where rj {= t — r) denotes 
the retarded time, rjo is an arbitrary constant, and H{x) the Heaviside func- 
tion which is 1 for a; > and for a; < 0. B{ri) is a smooth and monotonically 
increasing function [Such a requirement comes from the consideration of the 
energy-momentum tensor (EMT). See the discussions given below.]. Clearly, 
when B{ri) = 0, the metric represents the Minkowski spacetime written in a 
cylindrical coordinate system. When -8(77) 7^ 0, it represents an exploding 
cylindrical null dust cloud accompanied by a cylindrical gravitational wave. 
To see this, let us first consider the corresponding EMT, which is now given 
by [11] T^jy = B'{vi)H{rj — riQ)/{2r)kfj_ku, where k^ is a null vector defined as 
kfj_ = ^li~^li- A prime denotes the ordinary differentiation with respect to the 
indicated argument. It is clear that this EMT represents a null dust cloud, 
consisting of massless particles and moving along the null hypersurfaces 
1] = const. The presence of a gravitational wave can be seen by calculating 
the Weyl tensor, C^^xa, which is thought of as representing the purely grav- 
itational field. In the present case, it turns out that the only non-vanishing 



component of it is ^o = -C^uXaL'^M'^L^M'' = -e^'^'i^B'{r])H{r] - r]o)/{2r), 
where L^" = e^^''^/^^^/^ _ S^)/^ and M^ = {6^ + ir-^6!^)/V2, which means 
that the spacetime is Petrov type N, a typical gravitational wave spacetime 
[12]. In the present case, the gravitational cylindrical wave moves outward 
along the surfaces rj = const. From the expressions of T^j, and \E'o, we can 
see that the outgoing gravitational wave and dust cloud exist only after the 
time t = rjo, since H{ri — rjo) is zero for t < rjQ. 

On the other hand, it was shown in [13] that for a static CS the spacetime 
outside it is described by a flat cylindrical metric with a wedge of the angular 
size 6lp = Situ being removed, where n is the linear mass density per unit 
length. Thus, the superposition of the above two spacetimes will be the met- 
ric (1) with dip being replaced by {l~4fi)dip. Obviously, the resulting metric 
represents a cylindrical null dust cloud plus a gravitational wave propagating 
outwards in a CS background. Due to their presence the angle defect now is 



Sip = 2-n 



1 - (1 - Aii)r/ / e-^'^'^^/^dr 







(2) 



which is always less than Svr/i. This happens because both, the dust cloud and 
the gravitational wave, carry away positive energy. The latter can be seen 
from Thome's C— energy [14], which here is given by C = [1 — (1— 4/i)^e^*^'')]/8. 
Thus, the rate of the energy loss is C,^ = —(1 — 4/i)^e^*^^'-'6'(r7)/8, which is 
non-positive. That is, the outgoing gravitational wave and dust cloud always 



carry away positive energy and make the angle defect smaller. 

Taking the above solution as valid only in the region r > tq > 0, where 
Tq is the coordinate radius of the CS, and then extending it into the CS, we 
finally obtain a solution which can be written in the form of Eq.(l) but with 
the metric coefficient g^,^ being replaced by 

g^^ = -w^ = -P{rfH{ro - r) - (1 - Afxfr^H{r - tq), (3) 

where f3{r) is a smooth function with the boundary conditions /9(ro) = (1 — 
4/i)ro, and /9'(ro) = (1 — 4/i), which assure that the hypersurface r = tq is 
a regular boundary surface. When B{ri) = 0, the above solution reduces to 
the general static CS solution studied by Linet in [13]. When B{ri) ^ 0, it 
represents a CS, which is static for t < rjo and starts to emit gravitational 
waves and massless particles at the moment t = rjQ. To show this explicitly, 
let us first consider the corresponding EMT. It is easy to find 

Tf,u = (^{gfiu + r^Ty + (f^(fu) + pk^K, (4) 

a = -e^^^^jS'lS-^HiTQ - r), p = B'H{r] - r]o)w,r /{2w) 

w,r /w = (3'(3-^H{ro - r) + r'^H{r - ro), (5) 

where r^ = {—grrY^'^5^^,(p^ = {—g^ipY^'^^'^- Thus, providing that (3"{r) < 0, 
the first term in r.h.s. of Eq.(4) represents a CS with support in the region 
< r < tq. The regularity conditions at the axis r = require 

p{r) ^ r[l - 0{r^)], as r ^ 0. (6) 



To have the circumference monotonically increasing as r increases, we further 
require /3'(r) > (cf. Linet in [13]). The second term represents a null dust 
cloud moving outward from the axis. Since p is proportional to H{ri — rjo), 
we can see that this cloud comes to exist only after the moment t = rjQ. Note 
that the energy density of it is continuous across the hypersurface r = tq. On 
the other hand, from the equations T^x ''^ = we find T^^' = —{pk^kx)'^ = 
^ah'{vi)k^, where Tff^ is the EMT for the CS given by the first part in r.h.s. 
of Eq.(4), and a is given by Eq.(5). The above expression shows clearly that 
the CS is indeed the source of the emission. Outside the CS, T^^' = 0, 
and the dust cloud propagates to the space-like infinity with its amplitude 
decreasing like r"^. At a moment, say, t, the energy per unit length inside 
the radius r, is 

C = {l - [(3"H{r, - r) + (1 - ApfH{r - ro)]e^('')} /8, (7) 

from which we have that C,^ is always non-positive. That is, the CS is 
losing its energy through radiation. The maximum amount of the energy 
loss is Crnax = [1 ~ (1 ~ 4yu)^]/8 (~ /i wheu /i ^ 1). This gives a limit on 
the function h{ri) < — 2/n(l — Ap). When it is greater than this limit, the 
outgoing wave and cloud will carry away energy more than the CS has. Then, 
the CS will effectively have negative mass [cf. Eq.(7)]. In the following we 
shall not consider this possibility. 



7 



On the other hand, introducing a null tetrad frame by 

N, = {guf'kj^, M, = [i-g..)'^'5; - ti-g^^Y/'5^^]/V2, (8) 

where /c^ is defined as in (4), we find that the non- vanishing components of 
the Weyl tensor in this frame are 

^2 = -C^^xaL'^N'^iL^N'' - M^M")/2 = %% (9) 

where 

^^^ = ^- = -3^^^ = e^('')/57(4/5)iJ(ro - r), 

^'q" = -e''^'^^b'{r])w,r/{2w). (10) 

The reason for calculating the Weyl tensor in the null frame (8) is that its 
components in this frame have direct physical meaning, for example, \E'o 
represents the gravitational wave component propagating along the null hy- 
persurfaces defined by A*"^ ( or equivalently k^), and \1'4 the component prop- 
agating along the null hypersurfaces defined by L^ [12]. Equation (9) shows 
that ^0 contains two terms, the first is finite and has support only inside 
the CS, and so do the ^2 and \1'4 terms. We interpret them as the gravi- 
tational field of the CS. The second term of ^0 is proportional to b'{ri) and 



has support in the whole spacetime. Moreover, it is continuous across the 
hypersurface r = tq. Outside the CS, this term has a typical form of a gravi- 
tational wave. So, we interpret it as representing the gravitational radiation 
of the CS. Therefore, we conclude that the solution given by (1) and (3) in- 
deed represents the gravitational and massless particle emission of a CS. Note 
that although the CS is emitting gravitational waves and massless particles, 
Eq. (5) shows that the energy density of it still increases as the time develops 
[recall that b{ri) is a monotonically increasing function]. This is because, as 
the CS is radiating, it also undergoes a collapse. The proper radius of it now 
is given by Jq" e~^^^^^'^dr, which decreases as the time passes. However, since 
b{ri) is a function with a limit, it never collapses into a singular line. 

From Eqs.(l) and (3) we can show [10] that the null vector /c^ defines a 
null geodesic congruence, and that the quantity Qk = —k^;\= ~e''^^^w,r /w, 
represents the rate of the contraction of the congruence. Considering Eqs. 
(5) and (6) we find that Qk -^ — oo as r -^ 0, which indicates that when the 
backreaction of the radiation is taken into account, the axis might become 
singular. To show that this is indeed the case, let us turn to consider the 
Kretschmann scalar, which is found to take the form 

n = i?^,Aai?^'^" = 8[8(^;)^)2 + a^]/3 + 16^2'^r + 4e'(''V^, (H) 

where p, ^q*, \E'4, and '^q' are given, respectively, by Eqs. (5), (9) and (10). 



Equation (11) contains three terms, each of which has the following physical 
meaning: The first term represents the contribution of the CS, which is finite. 
The second represents the interaction of the gravitational field component \I'2* 
of the CS with the outgoing wave \E'o^, while the last represents the interaction 
between the matter field of the CS and the outgoing null dust cloud. From 
Eqs.(6), (9) and (10) we can see that the last two terms become unbounded 
as r — > 0. That is, because of the backreaction of both gravitational and 
massless particle radiation of the CS, the spacetime becomes singular at the 
axis. This singularity is a scalar one and always formed, as long as the 
radiation is not zero, characterizing by the function b{ri). A particular case 
is /3(r) = asm{6r), where a and 6 are two constants and determined by the 
boundary conditions at r = tq. When b{ri) = 0, the solution reduces to the 
one studied by Hiscock in [13]. By properly taking the limit tq — ^ 0, Hiscock 
obtained the solution of an infinitely thin CS with support only on the axis. 
Equation (10) shows that after taking this limit the \E'^'s are still different 
from zero and have support on the axis, too. This observation is important, 
since it indicates that our above results concerning the singularity behavior 
also hold for the case of an infinitely thin CS radiating gravitational waves, a 
case often considered in the numerical stimulations of long CS networks [3]. 
On the other hand, let us consider the timelike geodesies perpendicular 
to the CS. The first integration of the geodesic equations yields dt/dr = 
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(1 + e^) / 2 , dr / dr = (1 — e^)/2,dz/dT = 0, and dip/dr = 0, where r is the 
proper time. Perpendicular to the tangential vector AJ^qn (= dx'^/dr), we can 
construct other three linearly independent spacelike vectors A(^N(a = 1,2,3) 
by A(\-) = {dr/dT)5t + {dt / dT)5i^ , \%-. = 6^, and X'f^-. = {l/w)6^. One can 
show that these four unit vectors form a freely-falling frame. Computing 
the Riemann tensor in this frame, we find that it has only two independent 
components, one of them being given by 

i2M-^Afo)A[3)Afo)A(3) = [2e"'b'ir])w,r-{l ~ e''fw,rr]/{^w). 

The combination of the above equation with Eqs.(3) and (6) shows that as 
r — !► 0, this component diverges like r^^. In other words, the tidal forces felt 
by the freely-falling test particles become infinitely large at the axis. 

It might be argued that the formation of the singularity is due to the 
particular solution constructed above. However, the following considerations 
show that it is formed for all the radiating CS. To see this, let us consider 
the spacetime described by the metric [15] 

ds' = e2(^-^)(rft2 _ dr^) _ e'^'dz^ - e-'^^W^d^^ 

where K, U, and W are functions of t and r. The regularity conditions at the 
axis require 

\dv\' = -9^^ = e-^''W'^0{r'), as r ^ 0. (12) 
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Now, taking Eq.(4) as the source of the Einstein field equations, one can show 
that the null dust has only contributions to the metric coefficient K{t, r) (for 
a more general discussion on this point, see [10]). In fact, writing K{t,r) as 
K{t,r) = K'^'^{t) — b{ri)/2, the Einstein field equations can be divided into 
two groups, one reads 

p = b'{r]){W,r-W,t)/{2W), a = e^^'''''\W,u-W,rr)/W, (13) 

and the other consists of the equations for the functions K'^^, W, and U, which 
are given exactly by those for the case where b{ri) = and are provided in 
[15]. In addition to Eq.(12), we further assume that a is finite, so that when 
the radiation is switched off [b{ri) = 0] the spacetime is free of any singularity 
at the axis [16]. This assumption together with Eq.(12) yields 

W{t, r) -^ wir + wsr^ + 0(/), 

e^it'^) ^wi + air + 0{r^), as r ^ 0, (14) 

where wi,W3 and oi are smooth and bounded functions of t only. Note 
that the other Einstein equations may impose further restrictions on these 
functions, however, we ffnd that Eq.(14) is sufficient for our present purpose. 
Also, since we are mainly interested in the behavior of the spacetime near 
the axis, we shall focus our attention only in the inside region of the CS, 
and assume, without loss of generality, that a matching to outside is always 
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possible. With the above assumptions, we find that the Weyl tensor has 
only three non-zero components, \1'a {A = 0, 2, 4) in the null tetrad frame of 
Eq.(8), and can be written as ^^ = ^^ + ^^^, where 

n = -^/4 + fA, ^r = K = 0, 

^.. = _e2(^-^)6'(r^)[(f/„-f/„) - {W,,-W,r)/{2W)], (15) 

and j'a are functions of U and its upto second-order derivatives. From Eq.(14) 
we can see that they are finite at the axis. In sequel, the \E'^'s are finite, too. 
As before, we interpret these \l'^'s as representing the gravitational field of 
the CS, and ^q^ [which is proportional to b'{ri)] the gravitational radiation 
of it. In terms of these quantities, the corresponding Kretschmann scalar is 
given by 

n = 16[3(^^^)2 + %'%' + ay 6] + 16^2'^r + 4e2^^"^Vp, (16) 

where a, p, \I'^, and ^g^ are given by Eqs.(13) and (15). Similar to Eq.(ll), 
the above equation also contains three terms, each of which has the same 
physical interpretation. In particular, the last two terms represent the inter- 
action of the gravitational and matter fields of the CS with the corresponding 
radiation fields emitted by it, and behave like r^^ as r -^ [Note that here we 
still have T^^' = —{pk^k\)''^ = ^ab'{ri)k^.]. On the other hand, we also have 
Qk = -k^;x= -e^^^-^\W,t+W,r)/W ^ -r'^ -^ -cx), as r ^ 0. There- 
fore, it is concluded that the spacetime of a radiating CS is always singular at 
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the axis, because of the non-linear interaction of the gravitational and mass- 
less particle radiation fields emitted by the CS with the corresponding ones of 
the CS. 

The implication of our results to CS is obvious. In particular, it gives 
rise to another challenge to the structure formation scenario from CS, unless 
one assumes that the gravitational and particle radiation of ILCS is com- 
pletely suppressed so that the singularity will not be formed, or somehow 
considers the singularity as physically acceptable. In addition, "wiggle" CS 
have different properties from the "normal" ones, and the backreaction of 
the gravitational and particle radiation of them may not necessarily produce 
spacetime singularities, a problem which is now under investigation. 
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